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Abstract. Most calculations of air shower development have been based on the Bethe- 
Heitler cross sections for bremsstrahlung and pair production. However, for energetic 
^ ■ enough particles, a number of different external factors can reduce these cross sections 

OO \ drastically, slowing shower development and lengthening the showers. Four mechanisms 

0^ ■ that can suppress bremsstrahlung and pair production cross sections are discussed, and 

their effect on extremely high energy air showers considered. Besides lengthening the 
showers, these mechanisms greatly increase the importance of fluctuations in shower 
' development, and can increase the angular spreading of showers. 

Qh| Presented at the Workshop on 

Q ■ ^'Observing Giant Cosmic Ray Airshowers for > 10^^ eV Particles from Space" 

November 13-15, 1997, College Park, MD 

a 

INTRODUCTION 

^. 

, The electromagnetic portion of high energy air showers is governed by 

bremsstrahlung and pair production. Although the formulae for these process have 
been around for over 60 years, it is not well known that, in many situations, these 
formulae can be very wrong. The medium in which the bremsstrahlung or pair 
production occurs can drastically affect the cross sections. This contribution will 
discuss four different ways in which the medium can reduce the bremsstrahlung 
and pair creation cross sections. 

These suppression mechanisms can affect air showers by increasing the effective 
radiation length, lengthening the showers, and moving the position of the shower 
maximum deeper into the atmosphere. For ground based arrays like the proposed 
Auger Observatory, even a small change in the depth of shower max can affect the 
energy hitting the ground, especially for non-vertical showers where shower maxi- 
mum is considerably above ground level. Air fluorescence detectors like Flys Eye 



and the proposed OWL can measure the shower profile, and so their energy mea- 
surement would be less affected by unforeseen changes in the shower development. 

However, a change in the position of shower maximum can affect measurements of 
the composition of the highest energy cosmic rays. Moreover, these mechanisms 
wiU drastically reduce the number of particles in the early stages of the shower, 
changing the shower development profile. 



I BREMSSTRAHLUNG AND PAIR PRODUCTION 
SUPPRESSION MECHANISMS 

Suppression mechanisms for bremsstrahlung and pair production are possible 
because of the unusual kinematics in these processes. For ultrarelativistic particles, 
the momentum transfer between the radiating electron or converting photon and 
the target nucleus is very small, especially in the longitudinal direction [1]. For 
bremsstrahlung where E ^ m, 

q\\^Pe-p'e-k/c (1) 

^^{E/cr - {mcY-^{{E-k)/cY - (^^)^ - = (2) 

where and Pg are the electron momenta before and after the interaction re- 
spectively, k is the photon energy, m is the electron mass and 7 = E/m. For 
ultrarelativistic electrons, q\\ can be very small. For example, for a 1 EeV electron 
emitting a 100 PeV photon, q\\ = 10~^eV/c. Because q\\ is so small, by the uncer- 
tainty principle, it must take place over a long distance, known as the formation 
length: 

h 2hE{E - k) 

vo — — — TTi ■ 

g|| m'^cfk 

For the above example, Ifo is 200 meters; for a 1 PeV photon from the same elec- 
tron, Ifo rises to 20 km. This distance is the distance required for the electron and 
photon to separate to become distinct particles. It is also the path length over 
which the emission amplitude adds coherently to produce the emission probability. 
If something happens to the electron or nascent photon while it is traversing the 
formation zone, then the coherence can be disrupted, reducing the effective for- 
mation length and hence the emission probability. Even weak forces, acting over 
a long formation length, can be strong enough to destroy the coherence required 
for emission. The mechanisms discussed here work by disrupting the electron or 
photon, reducing the effective formation length. 



A Multiple Scattering (The LPM Effect) 



Multiple scattering can cause disruption by changing the electron trajectory. If, 
taken over Ifg, the electron multiple scatters by an angle larger than the typical 



bremsstrahlung emission angle I/7, then emission can be suppressed [1]. 

The reduction can be calculated by considering the effect multiple scattering 

has on gy; as the electron changes direction, it's forward velocity is reduced, and, 
with it, producing a change in q\\. This can be modelled by dividing the multiple 
scattering evenly between p^, and Pg. Then, 



q\\ = ^J{EcoseMS/2/cY - {mcY -^{{E- k) cos^m5/2/c)2 - (mc)^ - k/c (4) 

where 6ms/2 is the multiple scattering in half the formation length. Eg/ E^lf/2Xo, 

where Es = m\J ^lTxci = 21 MeV, and Xq is the radiation length. Scattering after 
the interaction is for electron energy E — k. This leads to a quadratic in //: 
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(5) 



If multiple scattering is small, this reduces to Eq. (3). When multiple scattering 
dominates 



2ncE{E-k)Xo . kELPM 
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where Elpm is a material dependent constant, given by Elpm = m^c' XQ/2hEl 
3.85 TeV/cm Xq. For lead, Elpm = 2.2 TeV, while for water Elpm = 139 TeV 
and for sea level air Elpm = 117 PeV. 

Since the formation length is the maximum distance over which the 
bremsstrahlung amplitude add coherently, the bremsstrahlung amplitude is pro- 
portional to the formation length, so the suppression factor is 



da/dk ^ If _ / kELPM 
~ duBH/dk ~ lfo~ \ E{E -k) ^ ^ 

and the dN/ dk \/k found by Bethe and Heitler changes to dN/ dk 1/ \fk. 

A similar effect occurs for pair production, where the produced electron and 
positron can multiple scatter. The two effects are closely related, as is shown in 
Fig. 1, and this relationship can be used to relate the bremsstrahlung and pair 
creation formation lengths and cross sections. For pair production 

2hE{k-E) ^ 



the corresponding suppression is 
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FIGURE 1. Schematic Representation of bremsstrahlung and pair conversion, showing the for- 
mation zone. 



This calculation has several limitations. The semi- classical approach may fail for 
k ^ E. The calculation neglects the statistical nature of multiple scattering, instead 
treating it deterministically. And, it neglects many niceties like the large non- 
Gaussian tails on Coulomb scattering and electron-electron inelastic scattering. 

Migdal developed a more sophisticated approach which avoided many of these 
problems [2]. He treated the multiple scattering as diffusion, calculating the av- 
erage radiation for each trajectory, and allowing for interference between nearby 
collisions. He found the cross section for bremsstrahlung is 



ia^ ^ ii^|(£){^,GM +2[1 + (1 - yms)}ZH. 

where G{s) and (f){s) are the solutions to differential equations, given by [3] 
(p{s) = 1 - exp - 6s[l + (3 - 7r)s] + sV(0-623 + 0.796s + 0.6585^) 

Tpis) = 1 - exp - 4s - 8sV(l + 3.96s + 4.97s^ - 0.05s^ + 7.5s^) 
G'(s) = 3V'(s) -2(/)(s). 

where 
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For /c<c£', s~l/ < ■j9ms >■ For s 3> 1, there is no suppression, while for 
s <^ 1, the suppression is large, ^{s) accounts for the increase in radiation length 
as photon emission drops. Migdals solution for s and ^(s) is recursive, because s 
depends on C(s). The recursion can be avoided by defining [3] 
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FIGURE 2. ydaLPu/dy for bremsstrahlung as a function of y for various electron energies in 
lead {Elpm=^-'^ TeV), showing how the spectral shape changes. Electrons of energies 10 GeV, 
100 GeV, 1 TcV, 10 TeV, 100 TeV, 1 PeV and 10 PeV are shown. The upper curve (soUd Une) is 
for E = 10 GeV, and the emission drops as energy rises. 

This is possible because ^ depends only logarithmically on s; a modified formulae 
for ^ depends only on s': 



scahng is recovered, albeit with a different coefficient. 

Fig. 2 shows how the LPM effect reduces yda/dy {y — k/E) for electrons in a 

lead target. The 10 GeV electron curve is very close to the Bethe-Heitler prediction; 
in the absence of suppression, this curve would hold for all electron energies. As 
the electron energy rises, emission drops. At the highest electron energies, photons 
with k <^ E are almost completely suppressed. For E'^/E^pm < k < I.SE'^/Elpm, 
the Migdal curve rises slightly above the unsuppressed; this is a consequence of 
either the approach or the approximations Migdal used. 

Fig. 3 shows how the pair production cross section is reduced. Compared with 
bremsstrahlung, pair production suppression sets in at higher energies. Symmetric 
pairs are suppressed the most; in the extremely high energy limit, one of the pro- 
duced electrons takes almost all of the photon energy. So, where the LPM effect is 
extremely strong, an electromagnetic shower becomes a succession of interactions 
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FIGURE 3. dai,PM/dy for pair production in lead, as a function of y for various photon energies, 
showing how the spectral shape changes. The solid line is for A; = 1 TeV, and the cross sections 
drop with energy; the other curves are for photons of energies 1 TeV, 10 TeV, 100 TeV, 1 PeV, 
10 PeV, 100 PeV and 1 EeV. 

where an electron emits a bremsstrahlung photon that takes almost all of the elec- 
trons energy, followed by a very asymmetric pair conversion, producing an electron 
or positron with almost all of the energy of the initial lepton. 

Two metrics for suppression effects in showers are the electron energy loss 
{dE/dx) and photon conversion cross section. Fig. 4 shows how these two markers 
change with energy in a lead target. As the particle energy rises above Elpm (2.2 
TeV for lead), the energy loss and conversion cross section fall, with bremsstrahlung 
affected at lower energies than pair conversion. 

The LPM effect was studied experimentally by the SLAC E-146 collaboration, 
who sent 8 and 25 GeV electrons through thin (0.07 % to 6% of Xq) targets of 
materials ranging from carbon to gold. Fig. 5 shows their data for carbon [4] , the 
material that is closest to air. Both LPM and dielectric suppression are needed to 
explain the data. 

Since E-146, there have been a number of additional calculations of 
bremsstrahlung with multiple scattering, using a variety of different approaches 
[5-8]. Several calculations have showed that, the non-Gaussian tail of large angle 
Coulomb scatters introduces additional term into the cross section; in the limit of 
large suppression [6-8] 
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This will reduce the magnitude of suppression in extreme conditions. Inelastic 
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FIGURE 4. The reduction in electron energy loss {J^{kdN/dk)dk) due to the LPM effect (solid 
line) and in the photon conversion cross section (dashed line) due to the LPM effect in a lead 
target. The LPM effect turns on at higher energies for photons. 

interactions, involving the atomic electrons should be treated separatedly from the 
elastic nuclear interactions, and a different potential should be used [6,8]. This may 
explain the poor agreement observed by E-146 between the data and the Migdal 
calculations for their light carbon and aluminum targets, especially for the 25 GeV 
data, around k = E"^ /Elpm- In hght of these results, calculations based on Migdal's 
formulae should be treated with caution, especially for light targets like air, where 
the inelastic form factor is important. 

B Dielectric Suppression 

Dielectric suppression occurs because photons produced in bremsstrahlung can 
interact with the atomic electrons in the medium by forward Compton scattering 
[9]. The photon wave function acquires a phase shift depending on the dielectric 
constant of the medium t{k) = 1 — {hoOpY/k"^ where oOp is the plasma frequency 
of the medium. With this substitution, the photon momentum {k/c) in Eq. (1) 
becomes ^fekjc^ and gy acquires an additional term [fiuipY /2ck. This leads to a 
reduced Z/, and a suppression factor 

^ = k^ + linupf 

This effect only applies for small y,y < Udie = i^p/m. For typical solids, ujp ~ 60 eV, 
so Udie ^ 10"*^. For y < ydie, the photon spectrum becomes da/dk ~ k, suppressing 



0.15 



0.1 



0.05 





0.15 



0.1 



0.05 



(a) 2% Xq Carbon 

25 GeV beam 



(b) 6% Xq Carbon 

25 GeV beam 




(c) 2% Xq Carbon 

8 GeV beam 



(d) 6% Xq Carbon 

8 GeV beam 



0.12 



0.1 





0.2 



10 

k (MeV) 



100 500 0.2 



10 
k (MeV) 



100 500 



FIGURE 5. Comparison of data from SLAC-E-146 with Monte Carlo predictions for 200 keV 
to 500 MeV photons from 8 and 25 GeV electrons passing through 2% and 6% Xq thick carbon 
targets. The cross sections are dN/d{logk)/Xo where N is the number of photons per energy bin 
per incident electron, for (a) 2% Xocarbon and (b) 6% Xocarbon targets in 25 GeV electron beams, 
while (c) shows the 2% Xocarbon and (d) the 6% Xocarbon target in the 8 GeV beam. Monte 
Carlo predictions are shown for LPM plus dielectric suppression (solid line), LPM suppression 
only (dotted line) and Bethe-Heitler (dashed line); all include transition radiation. 



the emission by {'^hLVpY / k"^ . The effect is also clealy apparent in the E-146 data in 
Fig. 5. 



C Suppression of Bremsstrahlung by Pair Creation, and 

vice- versa 

As Landau and Pomeranchuk pointed out, when //o becomes as long as Xq, then 
partially created photons can be pair create, destroying the coherence [1]. A simple 
calculation of this effect can be done by limiting to a maximum of IXq [10]. This 
suppression is visible (stronger than LPM and dielectric suppression) for 

nc 

For E > Ep, there is a 'window' in k where this mechanism applies: kp- — 
Xo(jjp/2hc < k < kp+ — 2ficE{E — k)/{XoEg). In this region, the photon spectrum 
is suppressed by k, and da/dk is constant. For k < kp_, dielectric suppression 
dominates, while for k > kp^. the LPM effect is dominant. Ep is 1.6 PeV in water 
or ice, and 42 PeV for air at sea level. The EpS are so similar because the varia- 
tion in Xq and cUp partially offset each other. In sea level air, the 'window' is 331 
MeV< k<3.0x lO'^'^E^ (with E in eV). 

There should be a similar effect where a partially produced electron or positron 
emits a bremsstrahlung photon. This possibility limits the coherence over Ijq. 
The strength of this suppression has not yet been calculated, but it should be 
comparable to that for pair production suppressing bremsstrahlung. 

This approach is overly simplistic. The 'hard cutoff' in // should be replaced by 
a probabilistic approach, where the photon interaction probability depends on the 
distance travelled. And, other suppression mechanisms will also be in effect, the 
radiation length will be longer than the naive IXq. An accurate calculation should 
include the interplay between the two reactions to arrive at an overall effective 
shower distance. Still, the above equation gives a reasonable estimate of when this 
effect needs to be considered. 

Fig. 6 summarizes these results, and shows that a 'simple' bremsstrahlung photon 
spectrum can have several different slopes for different photon energies. 

D Magnetic Suppression 

An external magnetic field can also suppress bremsstrahlung and pair creation. 
The bending caused by an external field acts just as does multiple scattering. 
The difference is that the magnetic field bending is quite deterministic, while 
the scattering angles are randomly distributed. The magnetic bending angle is 
6b = Ap/p = eBl f sin {9 b) / E where B is the magnetic field and 9b is the angle 
between the field and the electron trajectory. Emission is reduced if 6*^ > I/7; this 
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FIGURE 6. Schematic view of bremsstrahlung da/dk when several suppression mechanisms 
are present. Below kc-, dielectric suppression dominates. Between kc- and kc+, suppression 

due to pair creation is dominant. Between kc+ and Ulpm = E'^ /Elpm, the LPM effect is most 
important, while above k^pM the Bcthc-Hcitlcr spectrum is present. For E < Ep, pair creation 
suppression disappears and LPM suppression connects with dielectric suppression. 



happens for y < sui{9b) / B^.^ where i?c = vn?e' jU is the critical field strength 
[11]. 

The magnitude of the suppression can be found using an approach similar to 
that used by Landau and Pomeranchuk for multiple scattering. Because Qb ~ //, 
while Qms ~ the energy dependence is different. With the definition Eb — 
mBc/B sin {9b), the suppression factor is [12] 
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A magnetic field will also bend pair produced electrons and positrons. The same 
equation holds, except that k — E replaces E — k. For a particle trajectory per- 
pendicular to Earth's ~ 0.5 Gauss field, Eb ~ 45 EeV. So, for cosmic rays with 
energies above 10^° eV, magnetic suppression must be considered. 



E Emission Angles 

Besides reducing the interactions rates, these mechanisms also increase the an- 
gular spread of electromagnetic showers. When bremsstrahlung and pair cre- 
ation occur with a large opening angle, the formation length is shortened. For 
bremsstrahlung 

2nE{E - k) 



where dr^ is the angle between the photon and the electron trajectory. When 
is included in the calculations, the average emission angle rises from 9-y ~ I/7 

to 9^ ~ I/'S'T- If S is small enough, this spreading can dominate over multiple 
scattering in determining the angular spreading of the shower. If the scattering is 
taken over 1/2Xq (assuming 1 interaction per Xq, and half the particles are charged 
and hence subject to multiple scattering), this occurs for roughly S < 0.05. Most 
calculations of air and water showers do not include this angular spreading. It is 
likely to be most important for calculations for radio emission from showers in air 
and ice [13]. 

II SUPPRESSION IN SHOWERS 

Suppression mechanisms can affect showers by slowing their development, effec- 
tively increasing the radiation length. Because they all work by reducing If, the 
effects do not add; instead, q\\ must be summed, and the total suppression calcu- 
lated. Usually, multiple scattering is the most significant cause of suppression, so 
it will be the focus of this section. 

The effective increase in radiation length due to multiple scattering can be seen in 
Figs. 4; it shows how the area under the curve in Figs. 2 and 3 drops as the incident 
particle energy rises. For bremsstrahlung, energy loss is halved for electrons with 
E = 22Elpm, while the pair production cross section is halved for k = IOOElpm- 
Then, the radiation length in the first generation of the shower doubles; succeeding 
generations will show smaller effects. 

However, the effects go beyond simply lengthening showers. Because soft photon 
bremsstrahlung is the first reaction to be affected, the number of interactions will 
decrease more rapidly then the electron energy loss. So, the initial part of a shower 
will consist almost entirely of a few high energy particles, without an accompanying 
'fuzz' of lower energy particles. The initial shower development depends on a much 
smaller number of interactions. For example, a 25 GeV electron in sea level air will 
emit about 14 bremsstrahlung photons per Xq, while a 10^^ eV electron will emit 
only 3. It is worth noting that, these numbers are naturally finite because dielectric 
suppression eliminates the infrared divergence. Pair production is similar; the pairs 
become increasingly asymmetric. The higher energy lepton from a 25 GeV photon 
pair conversion takes an average of 75% of k; for a 10^^ eV photon, the average is 
more than 90%. 

Because of this, shower to shower fiuctuations become much larger. Misaki stud- 
ied shower development in lead and standard rock. For E ^ Elpm he found that 
the position of shower maximum was shifted to larger depths, and that the position 
of shower maximum varied greatly from shower to shower, and that the shower to 
shower variations overshadowed the average shower development [14]. 

In the limit E ^ E^pm, the initial part of an air shower becomes a succession 
of asymmetric pair production, where the higher energy of the pair loses most of 
it's energy to a single bremsstrahlung photon, re-starting the process. In short, the 



paradigm that successive generations of air showers have twice as many particles 
with half as much energy as the current generation fails completely. 

Ill AIR SHOWERS 

The composition of the highest energy cosmic rays is still a mystery. This article 
will consider two possibilities: protons (the most popular) and photons. In both 
cases, we will take 3 x 10^° eV as a standard energy. For incoming heavy ions, the 
effects arc greatly reduced because of the lower per particle energy, and neutrons 
can be treated as protons. Because these are toy models, the likely possibility of 
the photon pair converting in the earth's magnetic field will be neglected. 

Most current works have considered proton initiated showers. There has been 
disagreement as to whether the LPM effect is important in air showers. Capdevielle 
and Atallah found that it had a large effect on 10^^ eV and 10^*^ eV showers [15]. 
However, Kalmykov, Ostapchenko and Pavlov found a much smaller effect; for a 
10^° eV incident proton, the number of electrons at shower maximum decreased 
by 5%, while the position of shower maximum shifted downward by 15 ± 2 g/cvo? 
[16]. This shift is less than the error on a typical measurement of shower maximum 
in a single shower, but it can have a significant effect on composition studies. Of 
course, as experiments probe higher energies, suppression becomes more important. 
None of these authors considered the effects of fluctuations or other suppression 
mechanisms. 

Air showers studies are complicated by the fact that density, and hence Elpm and 
Ep depend on altitTidc and temperature. Ignoring temperature changes, pressure 
decreases exponentially with altitude, with scale height 8.7 km. For showers, it is 
convenient to work with column depth measured in g/cvo?. In an isothermal model, 
then Elpm = 117 PeV(Ao/A), where A is the column depth and Aq is ground 
level, 1030 g/cm^. Temperature will modify the relationship; with a temperature 
correction this Elpm is 2.25 EeV at 36 g/cw? (1 Xq) depth, and 1 EeV at 90 g/cm^ 
(1 hadronic interaction length. A). Neglecting temperature, ydie — 1-3 x 1Q~^A/Aq. 
For pair creation suppression, Ep — 42 PeVy' {Aq/A). The corresponding photon 
window is 331 MeV < A; < 3.0 x 1{)~'^'^E'^{Aq/A). The two window 'edges' have 
different A dependencies because all three mechanisms have a different dependence 
on Up and Xq. 

Incoming photons react by pair production, while protons interact hadronically. 
A central hadronic collision will produce a shower of several hundred pions; the 
neutral pions will decay to photons. The highest energy 7r° will have a rapidity 
near to the incoming proton, and their decay photons will have energies around 
2 x lO^^eV. Many diffractive processes, such as A production can produce photons 
with similar energies. Overall, photons from central interaction will have an average 
energy of about 2 x 10^^ eV. 

Although a complete Monte Carlo simulation is required to understand the effects 
of suppression in air showers, simple calculations can provide some indications 



where it matters, and should differentiate between the results of Capdevielle and 
Atallah and Kalmykov and collaborators. Because Elpm decreases with depth, 

in concert with the average particle energy, suppression mechanisms can actually 
become stronger as one moves deeper in the atmosphere. The solid curve in Fig. 7 
shows Elpm as a function of altitude. The dashed curve shows the average particle 
energy, for an idealized Bethe-Heitler shower from a 3 x 10^*^ eV photon. In each 
successive radiation length, there are twice as many particles with half the energy. 
The curve with the short dashes shows a similar cascade, from a 2 x 10^^ eV photon 
starting at lA. The electromagnetic interactions at a given depth are determined 
by the ratio of the two curves, which gives E/Elpm- For the photon case, the 
maximum suppression occurs around 75 g/cm^, where E ~ 80£^lpm- Electron 
dE/dx is reduced about 80%, and pair production cross section is down by 60%; 
the radiation length has more than doubled. For the hadronic case, the effect is 
smaller, and, of course, these high energy photons are only a small portion of the 
total shower. On the other hand, hadronic interactions are only partially inelastic, 
and the proton may carry a significant fraction of it's momentum deeper in the 
atmosphere, where suppression is larger. Because of the large variations in energy 
deposition depths, it is difficult to give more quantitative estimates. 

This model underestimates the effect of suppression, because, with suppression 
shower development is slower than the idealized model, further increasing the 
amount of suppression in the next stage. However, it does show that suppres- 
sion is important in photon shower, and in at least parts of proton initiate showers. 
The effect is clearly smaller than that predicted by Capdevielle and Atallah, but is 
consistent with Kalmykov and collaborators. 

Beyond the affect on average showers, fluctuations must be considered, because 
the cosmic ray energy spectrum falls as dN/dE ^ l/E^, it is important to under- 
stand the tails of the energy resolution distribution; without accurate simulations, 
showers whose energy is overestimated can easily skew the measured spectrum. 

Fluctuations can affect both ground based arrays as well as air fluorescence de- 
tectors that optically measure the shower development. For ground based arrays, 
suppression can change the relative position of shower maximum and the detector 
array. Although a sea level detector is near shower maximum for 10^° eV vertically 
incident proton showers, for non-vertical showers, where the detector is deeper, it 
may be significantly behind shower maximum. The angular spreading discussed in 
Subsection IE must also be considered. 

Air fluorescence detectors will observe far fewer particles in the early stages of 
the shower than current simulations indicate. LPM, dielectric and pair creation 
suppression must all be considered. Although LPM suppression affects the widest 
range of energies, the other mechanisms produce a larger reduction in cross section 
where they operate. For E > Ep, for example, dielectric suppression will reduce the 
number of bremsstrahlung photons with k < 331 MeV by two orders of magnitude. 
Pair creation suppression is more energy dependent, but, for example, at 200 g/cvo? 
depth, emission of 10^° eV photons from 10^^ eV electrons will be reduced by a factor 
of 10. Monte Carlo simulations are required to find the actual shower profile, but 
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FIGURE 7. Elpm (solid line), average particle energy for a 3 x 10^° eV photon shower (dashed 
line) and average particle energy for a 2x lO^^eV photon created at lA (short dashes). At altitudes 
where the average particle energy is larger than E^pm, suppression is important; with the ratio 
of the two energies determining the degree of suppression. 

the initial stages of the shower will be much less visible than current simulations 
indicate; neither the profiles produced by Capdevielle and Atallah nor by Kalmykov 
will be accurate. 

Even relatively late in the shower, there will be a reduction in the number of 
low energy particles. For example, where the average particle energy is ~ 10^'^ cV, 
LPM suppression will reduce the number of photons below ~500 MeV. So, there 
can be measurable effects even relatively close to shower maximum, if a particle 
counting detector is too far above shower maximum, then it may underestimate 
the size of the shower. 

Another way to visualize how suppression works is to consider the probability 
of photons penetrating deep into the atmosphere. Because, for a given photon 
energy, suppression increases with depth, high energy photons have a non-negligible 
probability of penetrating deep into the atmosphere. This can drastically change 
the development of photons showers. For proton showers, it can create small, dense 
subshowers within the main shower. Fig. 8 shows the interaction probability as 
a function of depth for a 3 x 10^° eV vertically incident photon. The sohd and 
dashed curves show the LPM and Bethe-Heitler cases respectively. With Bcthe- 
Heitler interaction probabilities, essentially all of the photons have interacted by 
2>Xq in depth, while with LPM interaction probabilities, the photon has a 7% chance 
of surviving to IOXq. 

For protons, the effect is smaller. The other curves in Fig. 8 are for a 2 x 10^^ 
eV photon, produced by a hadronic interaction at lA in depth. With Bethe-Heitler 
(short dashes), almost all have interacted by 6 Xq, while with LPM (dots) cross 
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FIGURE 8. Interaction probability for photons from cosmic rays, as a function of depth in 
the atmosphere. The soUd and dashed line show the interaction probability for a 3 x lO^^eV 
photon incident on the top of the atmosphere for LPM and Bethe-Heitler interaction probabilities 
respectively. The short dashed and dotted lines show the interaction probability for a 2 x 10^® eV 
photon produced in an interaction at a depth of lA, also for LPM and Bethe-Heitler interactions 
respectively. 

sections, it has a 1.2% chance of surviving to IOXq. The effect on the average 
shower is small, but, because of the steeply falling spectrum, it is necessary to 
consider even relatively atypical showers. 

Not considered here is magnetic suppression, which can dominate over other 
mechanisms for interactions in the upper 3 g/cm^ of the atmosphere. 

IV CONCLUSIONS AND ACKNOWLEDGEMENTS 

Multiple scattering, the dielectric of the medium, bremsstrahlung and pair 
production themselves and external magnetic fields can suppress bremsstrahlung 
and pair production, with bremsstrahlung of low energy photons the most sub- 
ject to suppression. Suppression also broadens the angular spread of emitted 
bremsstrahlung photons and produced pairs. 

For energetic enough particles, these mechanisms reduce electron dE/dx and 
photon pair production cross sections. When this happens, electromagnetic showers 
are lengthened, and shower to shower fluctuations become much larger. When 
suppression gets very large, then the shower angular development can be dominated 
by interactions, rather than multiple scattering. 

The importance of suppression in air showers depends on the incident particle 
type and energy. If the incoming particles are photons, then at 3 x 10^° eV, the 



effects are large. For 3 x 10^° protons, the effects on the average shower is smaller, 
but the effects of fluctuations may be considerable, particularly when particular 
detector are included. As experiments probe to higher energies, of course, the 
effects will grow larger. 

The calculations presented here show the need for a complete Monte Carlo, in- 
cluding all relevant suppression mechanisms, in order to correctly determine the 
shower profile for showers with energies much above 10^° eV. 

I would like to thank my E-146 collaborators for many useful discussions. This 
work was supported by the USDOE under contract number DE-AC-03-76SF00098. 
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